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HALFTONING METHOD AND SYSTEM 



BACKGROUND 

[0001] A digital image can be represented as a matrix of picture elements 
commonly referred to as pixels. Each pixel has a corresponding intensity value. For 
an eight bit "continuo us tone" gray scale image, that value can range from zero to 
two hundred fifty-five with zero representing white, two hundred fifty-five 
representing black, and the intermediate values representing varying levels of gray. 
[0002] The digital image must be processed before it can be printed. A 
conventional black and white laser printer forms the image on a sheet of paper or 
other media by placing clusters of dots at selected locations on the sheet. Because 
printers cannot print continuous tones - whether it's the many shades of gray in a 
grayscale image or the millions of colors in a color photograph - digital images must 
be converted to halftones. This process is called halftoning or dithering. 
[0003] Traditional halftoning methods change either the size of printed dots or the 
relative density of dots on the page. These two approaches are analogous to 
amplitude modulation (AM) or frequency modulation (FM) used in communications. 
For example, black and white continuous tone photographs can contain millions of 
shades of gray. When printed these shades of gray are converted to a pattern of 
black dots that simulates the continuous tones of the original image. Simply put, 
lighter shades of gray are composed of fewer or smaller black dots spaced further 
apart. Darker shades of gray contain more or larger black dots closely spaced. 
[0004] With AM halftoning, the density of dot clusters, defined as the number of 
clusters of dots per unit area, is fixed. Tone rendition is achieved by varying the size 
of each dot. The most commonly used AM halftoning algorithm is clustered dot 
screening. Cluster dot screening has the advantages of low computational load, 
stable dot formation, and good resistance to such printer artifacts as dot gain and 
banding. Thus, it is widely used in laser printers where a single isolated dot may not 
develop stably. However, there are several drawbacks to AM halftoning. First, there 
is the well known tradeoff between the number of gray levels that can be rendered 
and the size of a screen period. Second, the regular dot placement in highlight 
regions may be very visible, and can cause Moire artifacts when periodic patterns in 
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the image are similar to the clustered dot frequency. This can be problematic when 
rendering images scanned from printed material. 

[0005] With FM halftoning, the dot size is fixed and the tone rendition is achieved 
by varying the dot density. Commonly used FM halftoning algorithms include 
dispersed-dot screening, error diffusion, and search-based halftone methods such as 
direct binary search (DBS). Error diffusion is perhaps the most popular FM 
halftoning algorithm. Though it requires more computation as compared to 
screening, it is still very efficient. In general, FM halftoning achieves higher spatial 
resolution than AM halftoning and is free of Moire artifacts. However, it can lack 
the print stability required for electro-photographic printing. 

[0006] The characteristics of AM and FM methods are complementary to each 
other. However, traditional halftoning methods take an either/or approach. What is 
needed is a printing method and system that can take advantage of both methods 
while avoiding their drawbacks. 

DESCRIPTION OF THE DRAWINGS 

[0007] Fig. 1 is a schematic representation of a computing environment in which 

embodiments of the present invention may be incorporated. 

[0008] Fig. 2 is a block diagram illustrating components of an image forming 

device in which embodiments of the present invention may be incorporated. 

[0009] Fig. 3 is a block diagram illustrating control logic components in which 

embodiments of the present invention may be incorporated. 

[0010] Figs. 4-5 are block diagrams illustrating a halftone module and look-up 

tables according to embodiments of the present invention. 

[0011] Fig. 6 is a block diagram illustrating the logical components of a halftone 
module in action according to an embodiment of the present invention. 
[0012] Figs. 7 and 8 help illustrate a traditional error diffusion process. 
[0013] Fig. 9 is a flow diagram illustrating steps taken to generate a halftone 
according to an embodiment of the present invention. 

[0014] • Fig. 10 is a flow diagram illustrating steps taken to produce and print a 
halftone according to an embodiment of the present invention. 
[0015] Fig. 1 1 is a flow diagram further illustrating steps taken to calculate a dot 
placement indicator according to an embodiment of the present invention. 
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DETAILED DESCRIPTION 

[0016] Various embodiments of the invented halftoning method modulate both dot 
size (AM modulation) and dot density (FM modulation). In addition to modulating 
the dot size and density, the halftoning method can also explicitly control the size of 
dot clusters at different intensity levels. Experiments demonstrate that "clustered" 
AM/FM halftoning can achieve high spatial resolution, smooth halftone textures, 
good printing stability, and Moire resistance. 

[0017] The terms digital image and halftone are used throughout the following 
description. A digital image is an electronic matrix of pixels selected and arranged to 
reveal any combination of text and/or graphics. Each pixel has an intensity value. 
For example, in an eight bit gray scale image, the intensity value for each pixel can 
range from zero to two hundred fifty five. To the human eye, the difference 
between one shade of gray and the next is not perceivable, so an eight-bit gray 
scale image can be though of as a continuous tone image. To be printed, a digital 
image must first go though a halftoning (dithering) process. The processed digital 
image is referred to as a halftone. The printed digital image can also be referred to 
as a halftone 

[0018] The description that follows is broken into sections. The first section, 
labeled "Environment, " describes an exemplary computing environment in which 
embodiments of the present invention may be implemented. The second section, 
labeled "Components, " describes exemplary logical and physical elements used to 
implement various embodiments of the present invention. The next section, labeled 
"Operation," describe s exemplary steps taken to practice various embodiments of 
the present invention. The last section, labeled "Look-Up Tables," describes a 
method for designing look-up tables that can be used by embodiments of the 
invented halftoning method and system. 

[0019] Environment: Fig. 1 illustrates a printing environment 10 in which it would 
be advantageous to implement embodiments of the present invention. Environment 
10 includes computers 12 and 14, image forming device 16, and computer readable 
media 17. Computers 12 and 14 represent generally any computing device such as 
a desktop computer, laptop computer, or PDA (Personal Digital Assistant) that may 
send a print job to image forming device 16. Image forming device 16 represents 
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generally any combination of hardware and programming capable of processing a 
print job and producing an image such as text and graphics on a media sheet such 
as paper or transparency. Examples include printers, copiers, and facsimile 
machines as well as multi-function devices that provide all three functions. It is 
noted that where an image forming device functions as a copier, it can create a 
print job without the need for computers 12 and 14. 

[0020] Computer readable media 17 represents generally any media that can be 
used by or in connection with an instruction execution system such as a 
computer/processor based system or an ASIC (Application Specific Integrated 
Circuit) or other system that can fetch or obtain the logic from computer-readable 
media 17 and execute the instructions contained therein. While shown as a 
diskette, computer-readable media 17 can be any media that can contain, store, or 
maintain programs and data for use by or in connection with the instruction 
execution system. Computer readable media 17 can comprise any one of many 
physical media such as, for example, electronic, magnetic, optical, electromagnetic, 
infrared, or semiconductor media. More specific examples of suitable computer- 
readable media 17 include, but are not limited to, a portable magnetic computer 
diskette such as floppy diskette or hard drives, a random access memory (RAM), a 
read-only memory (ROM), an erasable programmable read-only memory, or a 
portable compact disc. 

[0021] Computers 12 and 14 and image forming device 16 are interconnected by 
link 18. Link 18 represents generally a cable, wireless, or remote connection via a 
telecommunication link, an infrared link, a radio frequency link, and/or any other 
connector or system that provides electronic communication between devices 1 2- 
16. Link 18 may represent an intranet, the Internet, or a combination of both. 
[0022] The term print job refers to a series of instructions directing image forming 
device 16 to produce images on one or more media sheets. The instructions may 
include directions to form text, graphics, or a combination of both. A print job may 
be generated by computer 12 or 14 or directly by image forming device 16 where 
image forming device 16 functions as a copier. The instructions may also include 
finishing directions such as directions to print multiple collated copies. 
[0023] Components: The logical components of one embodiment of the invention 
will now be described with reference to the block diagrams of Figs. 2-4. Beginning 
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with Fig. 2, image forming device 16 includes, among other components not 
shown, print engine 24, memory 26, and control logic 28. 

[0024] Print engine 24 represents hardware and programming capable of utilizing 
print job data to print images on media sheets. Where, for example, image forming 
device 16 is a laser printer, print engine 24 includes an optical scanner, a photo 
conducting drum, toner and a fuser. The optical scanner modulates a laser beam 
across the drum according to halftone print code generated by control logic 28. 
The scanned portions of the drum attract toner. Toner is transferred from the drum 
to a media sheet forming the desired image. The fuser permanently adheres the 
toner to the media sheet. 

[0025] Memory 26 represents generally any memory to which electronic data can 
be written and retrieved. For example, memory 26 may be random access memory, 
flash memory, a hard disk, or any combination thereof. Control logic 28 represents 
generally any program capable of processing a print job and sending halftone print 
code to print engine 24. As used here, halftone print code refers to electronic 
instructions capable of being processed by print engine 24 to produce a replica of a 
digital image. 

[0026] When received by image forming device 16, a print job is typically in PDL 
(Page Description Language) format. The PDL format describes the layout and 
contents of a printed page or pages. Modern versions of the PDL format are object- 
oriented, meaning that they describe a page in terms of geometrical objects such as 
lines, arcs, and circles. Control logic 28 is responsible for converting a print job into 
halftone print code that can be used by print engine 24. Processing includes 
rendering the print job into a halftone. 

[0027] Referring now to Fig. 3, control logic 28 includes rasterizer 30, compressor 
32, halftone module 34, and manager 36. Rasterizer 30 represents generally any 
program capable of rendering a print job into one or more digital images with each 
digital image representing a page or portion of a page to be printed. A digital image 
that defines the location and properties of each pixel can be quite large and can 
quickly deplete memory 26. Compressor 32 represents any program capable of 
compressing a digital image into a more manageable size. Compressor 32 is also 
responsible for decompressing. 
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[0028] Halftone module 34 represents a program capable of converting a digital 
image to a halftone capable of being printed by print engine 24 - the halftone 
defining the properties of each dot to be printed by print engine 24. Manager 36 
represents any program capable of writing a compressed digital image to memory 
26 and retrieving compressed digital image from memory 26 so that it can be 
decompressed by compressor 32 and forwarded to halftone module 34. 
[0029] Fig. 4 is an exemplary diagram illustrating the logical components of 
halftone module 34. As shown, halftone module 34 includes placement control 38, 
and size control 40. A halftone can be thought of as a two dimensional matrix of 
cells organized into rows and columns. Each cell may or may not contain a dot. 
Placement control 38 represents a program capable determining which cells of a 
halftone will contain dots. In other words, placement control 38 is responsible for 
FM modulation. Placement control 38 is also responsible for controlling the size of 
dot clusters. For those halftone cells determined to contain dots, size control 40 
represents a program capable of specifying the size and placement of each dot 
within its given cell. This can be referred to as AM modulation. Examples of 
particular techniques used by placement control 38, and size control 40 are 
described below. 

[0030] Halftone module 34 also includes dot cluster look-up table (LUT) 42, dot 
density LUT 44, dot size LUT 46, and look-up table control 48. Each of LUTs 42- 
46 include a number of entries with each entry corresponding to a pixel intensity 
value. Each entry includes data that can be used when converting a digital image 
to a halftone. Each entry in cluster LUT 42 contains data representing a cluster 
factor corresponding to a given pixel intensity value or range of pixel intensity 
values. Each entry in dot density LUT 44 contains data representing a dot density 
factor corresponding to a given pixel intensity value or range of pixel intensity 
values. A cluster factor and a dot density factor are used by placement control 38 
to help determine whether a given cell in a halftone will contain a dot and the size 
of a particular dot cluster. Each entry in dot size LUT 46 contains data representing 
a dot size factor corresponding to a given pixel intensity value or range of intensity 
values. 

[0031] A size factor is used by size control 40 to specify the size of a dot 
contained in a given halftone cell. Look-up control 48 represents generally any 
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program capable of locating an entry in each LUT 42-46 that corresponds to an 
intensity value of a given pixel. Look-up control 48 is also responsible for retrieving 
data from each located entry. Referring to Fig. 5, it is noted that LUTs 42-46 might 
not be part of halftone module 34 and may instead be stored in memory 26 or 
another suitable location accessible to LUT control 48. 

[0032] It is noted that LUTs 42-46 may be replaced with a single combined LUT 
(not shown). Each entry in the combined LUT would correspond to a pixel intensity 
or range of pixel intensities and would contain a cluster factor, a dot density factor, 
and a dot size factor. 

[0033] Referring now to Fig. 6, imagine a digital image made up of a matrix of 
pixels arranged in columns and rows. Each pixel can be represented as P(m,n) 
where m and n represent the respective column and row containing the pixel. A 
pixel's intensity value can be identified as l(m,n). When pr ocessing a digital image 
to create a halftone, the pixels making up the digital image are examined 
sequentially. For each pixel, P(m,n), look-up control 48 locates an entry in dot 
cluster LUT 42 corresponding the intensity value l(m,n) for that pixel and retrieves a 
cluster factor, C(m,n), from that entry. In a similar manner, look-up control 48 
retrieves a dot density factor, D(m,n), from dot density LUT 44 and a dot size 
factor, S(m,n), from dot size LUT 46. 

[0034] It is noted that instead of acquiring factors C(m,n), D(m,n), and S(m,n) 
from LUTs 42-46, placement control might calculate C(m,n) and D(m,n) as a 
function of the intensity of pixel P(m,n). Similarly, size control 40 might calculate 
S(m,n) as a function of pixel intensity. 

[0035] Using the retrieved factors, C(m,n) and D(m,n), placement control 38 
determines whether a dot is to be placed in a halftone cell corresponding to the 
given pixel, P(m,n). This determination is represented by a dot placement indicator, 
B(m,n). B(m,n) can have a value of zero or one. A value of one indicates that a dot 
is to be placed. A zero value indicates that no dot is to be placed. Using the 
retrieved factor S(m,n) and the variable B(m,n), size control 40 determines the size 
of the dot to be placed. This output, referred to above as halftone print code, is 
represented by the variable HPC(m,n) and is sent to print engine 24. HPC(m,n), for 
example may be a pulse width modulation (PWM) code used by print engines of 
laser printers produced by the Hewlett-Packard Company. 
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[0036] Dot placement indicator, B(m,n), can be calculated using an error diffusion 
algorithm. Figs. 7 and 8 help illustrate a traditional approach to error diffusion while 
examples of particular techniques used by placement control 38 are described in the 
following section with reference to Figs. 9 and 10. 

[0037] Fig. 7 illustrates a five-by-three matrix of cells 50 with cell (m,n) located in 
its center. Cell (m,n) corresponds to a pixel P(m,n) of a given digital image. As 
noted above, the value of B(m,n) determines whether or not a dot will be placed in 
cell (m,n). B(m,n) is set to one if the adjusted intensity values (discussed below) of 
P(m,n) meets or exceeds a threshold value. If not, B(m,n) is set to zero. 
[0038] Imagine, for example, that P(m,n) is a pixel in an eight bit gray scale image. 
Its intensity value can range from zero to two hundred fifty-five. A threshold value 
may be set at one hundred twenty seven. A dot is placed in cell (m,n) only if the 
adjusted intensity value for P(m,n) is at least one hundred twenty-seven. On an 
eight bit scale, a cell with a dot has an intensity value of two-hundred fifty five. A 
cell without a dot has an intensity value of zero. Assume that the adjusted 
intensity value for P(m,n) is one hundred forty-five x exceeding the threshold value. 
Consequently, a dot will be placed, and cell (m,n) will have an intensity value of 
two hundred fifty-five - substantially larger that the intensity level of P(m,n). This 
difference in intensity values is referred to as an error, e(m,n). The error can be 
calculated by the equation: e(m,n) = pixel intensity minus dot intensity. In the 
example, e(m,n) = 145-255 = -110. 

[0039] This error can be diffused by shifting portions of the error to cells adjacent 
to cell (m,n). As seen in the example of Fig. 7, varying portions of the error are 
shifted to the following four cells (m+ 1 ,n), (m-1 ,n + 1), (m,n + 1), and (m+ 1 ,n + 1). 
The portion shifted to a given cell depends upon a weighting factor, W. The 
number of weighting factors equals the number of cells to which the error is 
shifted. In Fig. 7 there are four weighting factors - W(1), W(2), W(3), and W(4) - 
and the error e(m,n) is diffused as follows: 

• W(1) * e(m,n) is added to the intensity value for (m-h 1 ,n); 

• W(2) * e(m,n) is added to the intensity value for (m-1 ,n + 1 ); 

• W(3) * e(m,n) is added to the intensity value for (m,n + 1 ); and 

• W(4) * e(m,n) is added to the intensity value for (m + 1 ,n + 1 ). 
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The sum of the weighting factors should equal one. For example W(1) might equal 
7/16, W(2) might equal 1/16, W(3) might equal 5/16, and W(4) might equal 3/16 
so that W1 + W2 + W3+W4=1. 

[0040] As the error, e(m,n), is diffused, the intensity value for certain cells 
adjacent to (m,n) are adjusted. Referring now to Fig. 8, before determining whether 
a dot will be placed in cell (m,n), the intensity value for P(m,n) is adjusted by adding 
the weighted errors from cells (m-1 ,n), (m-1 ,n-1 ), (m,n-1 ), and (m + 1 ,n-1 ). The 
adjusted intensity value for P(m,n) can be represented by the variable lA(m,n) where 
| A (m,n) = l(m,n)+W1 *e(m-1,n) + W2*e(m + 1 ,n-1 ) + W3*e(m,n-1) + W4*e(m- 
1 ,n-1). It is this adjusted intensity value that is used to determine the value of 
B(m,n) and whether a dot is to be placed in cell (m,n). 

[0041] Operation: The operation of embodiments of the present invention will now 
be described with reference to Figs. 9-1 1 . Fig. 9 is an exemplary flow diagram 
illustrating steps taken to produce a halftone according to an embodiment of the 
present invention. Figs. 10 and 1 1 are exemplary flow diagrams that help illustrate 
steps taken to generate and print a halftone according to an embodiment of the 
present invention. 

[0042] Starting with Fig. 9, for each pixel of a digital image, a pixel intensity is 
obtained (step 52). For a gray scale digital image, intensity is synonymous with 
gray level. Dot density is modulated according to pixel intensity (step 54). Dot 
cluster size is controlled according to pixel intensity (step 56). Dot size is also 
modulated according to pixel intensity (step 58). With reference to the examples 
above, steps 54 and 56 can be accomplished by placement control 38 (see Figs. 4- 
6). Step 58 can be accomplished by size control 40 using output generated by 
placement control 38. 

[0043] Fig. 10 illustrates steps taken to produce and print a halftone according to 
an embodiment of the present invention. Initially, a print job is generated (step 60). 
Referring back to Fig. 1, this step can be accomplished using software running on 
computer 12 to generate and send the print job to image forming device 16. 
Alternatively, where image forming device 16 is a copier, step 60 may be 
accomplished using image forming device 16 alone. The print job is then rasterized 
creating one or more digital images (step 62). For each digital image, a variable X is 
set to equal the number of pixel columns in the digital image while a variable Y is 
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set to equal the number of pixel rows (step 64). Variables m and n are each set to 
equal one (step 66). The variable m represents the column of a particular pixel being 
analyzed while n represents the row. 

[0044] The intensity value, l(m,n) is acquired for pixel, P(m,n) (step 68) and cluster 
factor C(m,n), density factor D(m,n), and size factor J5(m,n) are retrieved (step 70). 
Referring back to Figs. 4 and 5, steps 60 and 62 can be performed by LUT 
controller 48. Dot placement indicator, B(m,n) is calculated as a function of C(m,n) 
and D(m,n) (step 72). Step 72, can be performed by placement control 38 which 
sets the value of B(m,n) to one if a dot is to be placed and to zero if not. Next 
halftone print code HPC(m,n) is generated as a function of B(m,n) and S(m,n) (step 
74). Steps 72 and 74 will be discussed in more detail below. 
[0045] Halftone print code, HPC(m,n) is sent to print engine 24 (step 76) - the 
code indicates whether or not a dot is to be placed and the size of the dot if it is to 
be placed. The variable m is increased by one (step 78). It is determined if the 
variable m exceeds X (step 80). If not, the process jumps back to step 68 so that 
the remainder of row n can be processed. Otherwise, the variable m is set to one 
and the variable n is increased by one (step 82). It is determined if the variable n 
exceeds Y (step 84). If not, the process jumps back to step 68 so that row n can 
be processed. Otherwise, all columns and rows have been processed. 
[0046] Fig. 1 1 expands on step 72 of Fig. 10. As noted above, dot placement 
indicator, B(m,n), can be calculated using an error diffusion algorithm. Fig. 1 1 helps 
illustrate an example of an error diffusion algorithm modified according to an 
embodiment of the present invention. A clustered dot screen is provided (step 86) 
- the dot screen like the digital image having pixels arranged in X columns and Y 
rows. A threshold value (to compare against P(m,n)) is calculated in step 88 
according to the equation: T(m,n) = To + C(m,n)*Ds(m,n). To represents a default 
threshold value, and Ds(m,n) represents a density factor for dot screen pixel (m,n). 
To , for example, may be a fixed value or it may be a function of D(m,n). 
[0047] Using T(m,n) as a threshold value, dot placement indicator B(m,n) is 
calculated according to an error diffusion algorithm (step 90). As discussed above, 
the value of B(m,n) (one or zero) determines whether a dot will be placed in cell 
(m,n). B(m,n) is set to one if the adjusted value of density factor D(m,n) equals or 
exceeds threshold T(m,n). Otherwise, B(m,n) is set to zero. 
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[0048] D(m,n) is adjusted by adding to it weighted errors resulting from the 
processing of prior pixels of the digital image. As an example, when processing 
pixel (1,1) no other pixels have been processed. C(1,1) and D(1,1) are acquired 
according to the intensity value for P(1,1). Because no prior pixels have been 
processed, D(1,1) is not adjusted. T(1,1) is calculated and it is determined if D(1,1) 
equals or exceeds T( 1 , 1 ) . If it does, B( 1 , 1 ) is set to one. Otherwise B( 1 , 1 ) is set to 
zero. An error, e(1 ,1 ), can then be calculated by the equation: e(1 ,1 ) = D(1,1)- 
B(1,1). That error, e(1,1) can then diffused to one or more other pixels yet to be 
processed. 

[0049] Referring back to the traditional error diffusion example discussed with 
reference to Fig. 8, D(m,n) might be adjusted by adding weighted errors resulting 
from the processing of four prior pixels of the digital image according to the 
following equation: D adjusted (m,n) = D(m,n) + W1 *e(m-1 ,n) + W2*e(m + 1 ,n-1 ) + 
W3*e(m,n-1 ) + W4*e(m-1 ,n-1 ). The values of each weighting factor may, for 
example, be fixed or dependent upon the value of D(m,n). 
[0050] With the value of B(m,n) set, an error, e(m,n), can be calculated and 
recorded (step 92). The recorded error (e(m,n) = Dadj U sted(m,n) - B(m,n)) can then 
be used when processing subsequent pixels. It is noted that where the choice of 
values for dot placement indicator, B(m,n), are zero and one, the value of a dot 
density factor, D(m,n), will be within a set of values ranging from zero to one. 
[0051] With dot placement indicator B(m,n) set, halftone print code HPC(m,n) can 
be calculated in step 74 (Fig. 10) according to the equation: PC(m,n) = S(m,n) * 
B(m,n). HPC(m,n) is sent to print engine 24 in step 76, and the next pixel is 
processed by repeating steps 68-74 of Fig. 10. 

[0052] Look-Up Tables: Using the process described above, a pixel's inten sity 
value, l(m,n), is used to acquire a cluster factor, C(m,n), a dot density factor, 
D(m,n), and a dot size factor, S(m,n). These factors, acquired from look-up tables 
(LUTs 42-46 shown in Figs. 4 and 5, are then used to generate a halftone print 
code HPC(m,n) for that pixel. These look-up tables, dot cluster LUT 42, dot density 
LUT 44, dot size LUT 46, should be designed so that a desired level of absorptance 
can be obtained for each possible pixel intensity value, l(m,n). The tables can also 
be designed to help optimize a variety of printing attributes including print quality 
and print stability. Therefore, it is possible to achieve improved quality rendering of 
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a desired absorptance by varying the various factors, C, D, and S, within the three 
LUTs 42-46. 

[0053] There are two important quantities that can be considered when designing 
LUTs 42-46. The first quantity is the output tone level, T(C,D,S), at cluster factor 
C, dot density D, and dot size S. Here we assume that T(C,D,S) is linear in 
absorptance, so that T(C,D,S) = 1 is a perfect black and T(C,D,S) = 0 is white. 
Intermediate values of T(C,D,S) vary linearly with luminance Y. The other quantity 
is print distortion, U(C,D,S) which specifies the print quality of intensity level 
U(C,D,S) rendered with the triple (C,D,S). The particular definition of U(C,D,S) 
depends on the specific printing application and the printing artifacts to be avoided 
such as banding and moire artifacts. An overriding objective is to select the triple 
(C,D,S) for each input intensity level that gives the best print quality subject to a 
constraint that the tone T(C,D,S) is correct. 

[0054] For each set of parameters, C and S, a tone curve can be measured as a 
function of the dot density D. This is done by printing a test pattern of gray 
patches in randomized spatial positions with values of D ranging from 0 to 1 and 
constant values of the parameter pairs (C,S). The tone of each patch is then 
measured to obtain the function T(C, . ,S). Let A(l) be the desired tone curve 
where I is an intensity level between 0 and 255. Then, for each parameter pair 
(C,S) the tone correction curve Rcs(l) is computed with the property that: 

T(C, Rcs(l), S) = AO) 
Notice that for each fixed dot size S and cluster factor C, T(C,D,S) is a monotone 
increasing function of the dot density D. Therefore, T(C,D,S) can be inverted with 
respect to D. We denote this inverse function by Tc,s" 1 (.). A dot density curve can 
be computed according to the equation: 

Rc,s(l) = Tcs^Ad)) 

This dot density curve is referred to as the tone compensation curve for dot size S 
and cluster factor C. In practice, some absorptance levels may not be achievable for 
certain dot size values S, even when the dot density is chosen to be its maximum 
value, Dmax. In this case A(l) exceeds T(C, Dmax, S), and Res is set to equal Dmax. 
[0055] After determining Rcs(l), the print distortion of rendering intensity level I 
with the triple (C, Rcs(l), S) can be observed. Let (C^S 0 *) where j and i each range 
in value from zero to N, represent the set of possible cluster and size factor values 
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from which to choose. For each pair (C (j) ,S (j) ) / a gray ramp is printed using the 
parameter triple (C U) , Rco>,s(j>, S U) ). This gray ramp is designed to achieve the desired 
tone curve A(l). For each intensity value, I, the values of j and i are selected that 
produced the best print quality at that intensity level. The indices j and i 
corresponding to this best quality output are designated as i(l) and j(l). These two 
indices are a function of intensity level I because, in general, the best choice of C 
and S will be different for each value of intensity level I. 

[0056] One approach to select the best values of C and S is to visually inspect the 
print samples and selecting the ramp with the best appearance. This subjective 
evaluation can be based on the following objective considerations. First, the printed 
patch should have good print stability. More specifically, this means it should be 
resistant to the banding artifacts and intensity level discontinuities. Second, the 
printed patch should be smooth with an absence of grain. Third, the halftoning 
algorithm at each gray level should be resistant to Moire artifacts. 
[0057] The following guidelines can be used in selecting the values of C and S that 
would achieve these objectives. In the highlight regions, the dominant artifacts are 
visible isolated dots. Therefore, it is generally desirable to use smaller dot size in 
these regions. Moreover, C = 0 is preferred for highlight regions in order to retain the 
"blue noise" d ispersed-dot distribution. For the midtone and shadow areas, print 
stability artifacts, such as banding, become more dominant. Thus, we prefer larger 
C values to generate larger dot clusters which are more stable for electro- 
photographic printing. However, if a C value is too large, it can compromise Moire 
resistance. Therefore, a tradeoff must be made between print stability and Moire 
resistance. To further enhance the stable formation of dot clusters, the following 
rule can be applied for dot size selection: if the cluster factor C>0, the dot size 
factor S should at its maximum. 

[0058] Once the values of |(l) and j(l) are selected, cubic spline interpolation is 
used to generate the dot size curve Si and the cluster curve Ci for intensity values 
ranging from zero to two hundred fifty-five. The dot density curve is computed as 
Di = Rci,si(l). It can be beneficial that the cluster, dot density, and dot size curves 
(Ci. Di, and Si) each vary smoothly in order to avoid contouring artifacts when 
rendering slowly varying image regions. 
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[0059] Conclusion: The diagrams of Figs. 2-6 show the architecture, functionality, 
and operation of various embodiments of the present invention. A number of the 
blocks are defined as programs. Each of those blocks may represent in whole or in 
part a module, segment, or portion of code that comprises one or more executable 
instructions to implement the specified logical function(s). Each block may 
represent a circuit or a number of interconnected circuits to implement the specified 
logical function(s). 

[0060] Also, the present invention can be embodied in any computer-readable 
media for use by or in connection with an instruction execution system such as a 
computer/processor based system or an ASIC (Application Specific Integrated 
Circuit) or other system that can fetch or obtain the logic from computer-readable 
media and execute the instructions contained therein. "Computer-readable media" 
can be any media that can contain, store, or maintain programs and data for use by 
or in connection with the instruction execution system. Computer readable media 
can comprise any one of many physical media such as, for example, electronic, 
magnetic, optical, electromagnetic, infrared, or semiconductor media. More specific 
examples of suitable computer-readable media include, but are not limited to, a 
portable magnetic computer diskette such as floppy diskette or hard drives, a 
random access memory (RAM), a read-only memory (ROM), an erasable 
programmable read-only memory, or a portable compact disc. 

[0061] Although the flow diagrams of Figs. 9, 10 and 1 1 show specific orders of 
execution, the orders of execution may differ from that which is depicted. For 
example, the order of execution of two or more blocks may be scrambled relative to 
the order shown. Also, two or more blocks shown in succession may be executed 
concurrently or with partial concurrence. All such variations are within the scope of 
the present invention. 

[0062] The present invention has been shown and described with reference to the 
foregoing exemplary embodiments. It is to be understood, however, that other 
forms, details, and embodiments may be made without departing from the spirit and 
scope of the invention that is defined in the following claims. 
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